Phonon superradiance and phonon laser effect in nanomagnets 
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We show that the theory of spin-phonon processes in paramagnetic solids must take into account 
the coherent generation of phonons by the magnetic centers. This effect should drastically enhance 
spin-phonon rates in nanoscale paramagnets and in crystals of molecular nanomagnets. 
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The problem of spin-lattice interactions is almost as 
old as quantum theory of solids pj. It is about comput- 
ing the amplitude of the transition between spin states 
of a magnetic atom due to the spontaneous emission of 
a phonon or due to, e.g., the Raman scattering of an 
existing thermal phonon. These processes are respon- 
sible for the width of the paramagnetic resonance and 
for other spin-lattice relaxation effects. In this Letter 
we shall demonstrate that all previous works on the sub- 
ject have missed the collective nature of the relaxation 
when the wavelength of the phonon is comparable to the 
distance between the magnetic atoms. General theoret- 
ical treatment of this effect has become possible due to 
the recently established mechanism of the dynamic spin- 
phonon coupling |2|, 



H s _ ph = —HS ■ fi, 



n(r) =-V x u(r), (1) 



where S is the spin of the atom, CI is the angular veloc- 
ity of the local rotation of the crystal due to a transverse 
phonon, and u is the phonon displacement field. Eq. (JTJ 
has the following origin. The crystal field (i.e., the mag- 
netic anisotropy) that determines spin states of an atom 
in a solid is defined in a local coordinate frame coupled 
to the crystal axes. In that coordinate frame the local 
rotation of the lattice is equivalent to the magnetic field, 
B = Ml/g(i B (with g being the gyromagnetic factor for 
the spin), that has Zeeman interaction with S. Since the 
local rotations of the crystal lattice by phonons are small, 
all matrix elements computed in the rotating frame and 
in the laboratory frame are practically identical. The 
beauty of this approach is that it has no unknown pa- 
rameters. For, e.g., spontaneous emission of a phonon, 
all parameters of the crystal field enter the interaction 
through the phonon frequency only. The latter is deter- 
mined by the distance between the spin levels and can be 
directly measured, leaving no fitting parameters for com- 
parison between theory and experiment. Another way 
to look at this interaction is through the conservation of 
the total angular momentum, spin + lattice. In order to 
conserve the angular momentum, the transition between 
the spin states of an atom produces a local elastic twist 
that propagates out as a transverse phonon. This Let- 
ter is based upon observation that an assembly of the 
magnetic atoms in the excited state must interact via 
transverse phonons very much the same as atoms inter- 



act via photons in the laser physics. The deep analogy 
between the two problems can be seen by noticing from 
Eq. that ii/2gfi B is equivalent to the vector poten- 
tial of the electromagnetic field. The huge (c/vt) 3 ratio 
of the phonon density of states to the photon density of 
states (with v t being the speed of the transverse sound) 
makes such a phonon laser much more powerful than the 
photon laser based upon magnetic dipoles. 

We consider two nearly degenerate states of spin S, 
described by the Hamiltonian for pseudospin 1/2, 



H 



1 . 



-hu>o ■ cr, Hujq = We z + Ae x , 



(2) 



where <x is the Pauli matrix, A is the splitting of the 
levels at the avoided crossing, and W = E-\ — E\ is the 
energy difference for the two degenerate states at A = 
(±1 denoting spin up and down). We label these states 
as j^-i) an d then E±\ = (V'iil -^s |^±i) an d A = 
2 (t/'-i Hs , where Hs is the crystal-field Hamilto- 
nian for the spin S. Choosing the new z axis for pseu- 
dospins in the direction of u>q, one obtains the eigenvalues 
and the eigenfunctions of H : E± = ±(l/2)\ / W 2 
and \ip±) — (C: 



i± w/Va 2- 



A 2 

^\i) 1 )±C±\i)_ 1 ))/ y/2 with C± = 
W 2 , whereas 

The projection of the spin-phonon interaction, Eq. ifTjl. 
onto the two-state basis is 

^ ^ «— » ^ 

H s ^ p h = ~^cr- g • (V x u) = --o- a g al3 (V x u)^ . (3) 

For tunneling between up and down states of spin S, 
\yj ±1 ) = \±S) and 



lx = e z SC+C- =e z S 
.1 



-e 2 S- (C 2 



^W 2 + A 2 ' 
C 2 ) = -e z S- 



Sy = 



w 



Vw 2 



(4) 



For many two-state particles coupled to phonons the 
effective many-body Hamiltonian is 

H = -fX^O'^IXy^' I v x H*i)]+Hph- (5) 



2 



where H p h is the Hamiltonian of non-interacting 
phonons. We use canonical quantization of phonons, 



2MN 



E 

kA 



ekAe' 



(«kA + aL kA ) 



(6) 



where M is the mass of the unit cell, N is the number of 
cells, e kA are unit polarization vectors, A = t, t, I denotes 
polarization, and uux — v\k is the phonon frequency. 
The spin-phonon Hamiltonian then becomes 



ph 



^5> ikr «o- r G kA (a kA 

kA i 



i -kA 



(7) 



where G kA is due to transverse phonons only, A = t,t, 

G kA = 



kX 

2MN 



g • [k x e kA ] , 



G kA = G 



kA- 



(8) 



In the Heisenberg representation for time-dependent 
operators, ih[dO(t)/dt] = [0(i),H], one obtains, with 
account of the commutation relations, [<Ji a ,<7jf}] — 
2ie a 1 a ij Sij and [a kA ,a k , A ,] = SwS\y, the following 
coupled equations: 



(Ti = (Ti x u> 



'G kx (a kx -al kx )} (9) 

OkA = - {i^k\ + 7/,-a) fl kA + ~ E e ~ 4k ' ri °V G kA ■ (10) 



kA 



Due to the linear interaction between cr.j and u, Eq. (0) 
conserves (Y2 i Gij 2 for uniform rotations. Consequently, 
(Ti in Eqs. @ and (|1(J|> can be treated as a semiclassical 
average over distances large compared to the spacing be- 
tween magnetic atoms but small compared to the phonon 
wavelength. We are interested in the interaction of such 
local polarization with a semiclassical phonon field. One 
can write the Heisenberg operators as 



er ,;+(£) = a l+ {t)e~ 



aux(t) = a kA (i)e~ 



(11) 

where &i+(t) and a kA (i) are slow variables. Keeping only 
resonant terms in Eqs. 10 and (|10|l . and neglecting the 
hybridization of phonon and pseudospin modes, one ob- 
tains after integration on time 



a kA 



4 uj k \ - cjo 



a i+ G- : 



kA, 



J 7fcA 



(12) 



where G kA= t = G kA x ± iG kA4) . The insertion of a kA and 



-kA 



into Eq. @ then gives the system of equations: 



7°J± 



2 aiz E r *. 



with 



27. 



w J 



■t2 



(13) 



(14) 



and 



N 

T 



E 

\=t,t 



V.X.X 



GkA.y) 



(15) 



For a small sample of size, L < Ao ~ Vt/ujQ, containing 
A/" S> 1 magnetic atoms or molecules (1 « ^ « N) 
that are initially polarized in one direction, one can use 
gjk-frt-rj) ^ i ^ Eq Eq S then preserve the 

length of the total polarization, Y^ , a j . The unit vector 
er = (l/J\f)(Y2j Cj) satisfies the Landau-Lifhitz equation, 



\<T X W 



a [cry. cr x uj \ 



(16) 



where a = 7VTi/(2cjo) is the dimensionless damping 
constant and Ti = TV -1 ^ k V£2ir5 (uJk — ^o) is the 
one-atom spin-phonon rate. Replacing ./V -1 ^ k ... by 
vq J d 3 k / (2ir) 3 . . . (vq being the unit cell volume) and 
writing V£ with the help of Eqs. |@J|, I©, and as 



v3 = 



8h Mv\ ul 



h 2 

1 _ _£ 

fc 2 



one obtains for the one-atom spin-phonon rate 0: 
1 ~~ 127rfi, Mv 2 oj 3 d ~ \2irh pv\ ' 



(17) 



(18) 



1 /3 

where ujd = Vt/v is the Debye frequency for the tran- 
verse phonons and p — M/vq is the mass density of the 
crystal. It is easy to check that Tsr = A/Ti is the rate of 
the relaxation of a z in Eq. I|ltj|) . that is, the longitudinal 
relaxation rate. The transverse relaxation rate for a± is 
Tsr/2. This is Dicke superradiance Q with the rate ex- 
ceeding the one-atom spin-phonon rate by a large factor 
N. Eq. i|TT)|) . among other situations, describes the case 
when all pseudospins are initially prepared in the excited 
state, a z — —1, and then relax collectively towards the 
ground state, u z = 1, by emitting coherent phonons of 
frequency wo- Because V£ of Eq. lfTT)l is proportional to 
sin 2 # = 1 — k 2 /k 2 , the angular distrubution of the emit- 
ted phonons is maximal in the direction perpendicular to 
the quantization axis z. 

For large samples, A < L, the phase of the emitted 
phonons is no longer constant throughout the sample, 
and one should also study inhomogeneous solutions of 
Eqs. (| 1 31) . For that purpose it is convenient to switch 
to the Fourier harmonics Oi± = A/" _1 5TJ q e±I q l " icr q± 
and <Ji Z — M^ 1 Y2q el ql " icr qz- When the paramagnetic 
body is a small part of a large solid matrix, the phonon 
modes are practically continuous in the k-space while q 
is quantized. The first of Eqs. ifl"^ then becomes 

CTq± = TiWoOq-t - \ E r ( C l") Cr [±(q-q")],z Cr q"± , (l 9 ) 



q" 



where 



r (q) = ^E^- q |2 



(uj k - ^o) 2 + l 2 k 



(20) 



3 



and F k = TV -1 J2i e * k ri — V 1 Jdre ikr . For a rectan- 
gular body of dimensions L a (a — x, y, z) it gives 



n 



sin(fc Q L Q /2) 
k a L a /2 



(21) 



whereas for a sphere of radius R one obtains — f(x) = 
3a; -3 (sin X — a: cos a;) with x — kR. 

Superradiance corresponds to the q = mode in Eq. 
ijTSfr. For 7 k — > it occurs at a rate 

Tsr - A^ / ——V£\F*\ 2 2K6{w k -u ) 



(22) 



where (. . .) denotes the average over the directions of ko 
that satisfy Wfc = ujq. For fc a L a <C 1 one can replace F^ 
by 1, and the second of Eqs. reduces to the Dickc 
rate, TgR = A/Ti. In the general case, the enhancement 
of the relaxation rate due to superradiance is given by 

Tsr n M Af 

~ = "fcT ' nM = V 



(23) 



where um is concentration of magnetic atoms and 
klV (vl\F ko f)/(Vl) 



T 



(24) 



depends on the wavelength and on the shape of the 
magnetic body. Since the maximal value of T is of 
order 1, the maximal enhancement is always given by 
the dimensionless ratio JiM/kg. For a sphere one ob- 
tains T = (Air/ 3) (k R) 3 f 2 (k Q R). The destructive inter- 
ference, T = 0, occurs at kgR = nn, where n is an inte- 
ger. At k R > 1, T(k R) ~ l/(k R). That is, r SR goes 
down to Ti only at koR ~ Um /kg, that can be very large. 
For a box of dimensions L a = b a L, the geometrical fac- 
tor T must be computed numerically. For a cube of size 
L, T(fco£) is qualitatively similar to that of the sphere. 
It almost vanishes at koL = 27rn, which corresponds to 
the destructive interference of phonons emitted perpen- 
dicular to the cube faces. The choice of L x ^ L y causes 
a more complicated pattern shown in Fig. 

An interesting observation is that the increase of the z 
dimension of the body beyond the wavelength of the ra- 
diation does not lead to the decrease of the superradiant 
rate. The latter is the consequence of the suppression 
of phonons emitted in that direction, Eq. 117J1 . For a 
sample of a prism shape with k$L x ~ k$L y -C 1 and 
arbitrary koL z , one obtains T = (k$L x ) (koL y ) T (koL z ), 
where T (u) = u at u = k$L z <C 1 and T (u) = 3ir/2 at 
u> 1, see Fig.^3. Another interesting case is a flat sam- 
ple with k^L x -C 1, by ~ b z ~ 1 and arbitrary u = kf)L. 
In this case one can write T = (koL x ) T (u, b y , b z ), where 
f = u 2 b y b z at u < 1 and f S 37r at u > 1, Fig. [TJd. 
The latter is a thin-film limit for which T depends nei- 
ther on b VtZ nor even on the shape of the film area. The 
dependence of T on the thickness of such yz-film is 

1 - cos (k L x ) 



T = 6tt- 



koL 



(25) 
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FIG. 1: a - T for rectangular geometry; b - T for the film 
and prism geometries. L a ~ b a L, a — x,y,z. 



In this geometry the phonons are emitted in the direction 
x perpendicular to the film. The emission rate can be 
as large as r S R, ma x ~ 13.7 (n M /kl) T x at k Q L x w 2.33, 
while the intensity of sound is proportional to the film 
surface. One should notice the difference from the case 
of the xy-film, Fig. for which the superradiant rate 
decreases with increasing the dimensions of the film sur- 
face beyond the phonon wavelength. 

Let us now consider the relaxation of the non-zero cr q 
harmonics in macroscopic samples for a particular case 
of the maximal inversion, a z = — 1. At k$L 3> 1, because 
of the factor |Fk-q| 2 in Eq. (|2Ufl. the phonon momenta 
k that contribute to Eq. H2U|) must be very close to q. 
Under the above conditions, Eq. I|19l) simplifies to 



1 , . 

cr q ± = T^oOqi + 2 r ( c l) cr q±' 
where T(q) is now the increment rate given by 

2 7 g 



r(q) = v 2 - 



\2 

u>o) 



(26) 



(27) 



The highest rate corresponds not to q = 0, as for the su- 
perradiance, but to q that provides the minimal value of 



4 



Lu q —iOQ determined by the previously neglected hybridiza- 
tion of the spin and phonon modes, (Aw) m j n ~ Vk . This 
is the laser mode that grows at a rate 



ko, max7fco 



v3 

ko , max 



ko , max 



S 2 A 2 m 
8h Mv 2 



(28) 



where V, 2 is the maximal value of V, 2 achieved for 

Ko , max Ko 

k perpendicular to the z axis. The maximal value of 
on 7 fco is TL, m ax = Hc ,max, and the enhancement of the 
spin-phonon relaxation is very large: 



L.max 



3^2 



~A~ 



1/2 



Mv 2 



1/2 



(29) 



At Vk ,max "C 7fc , Eq. (|2*8|) reduces to the well-known 
laser formula, Tl = V£ max /7fc - As in optical lasers, 
the increment rate Yl must compete with the dephas- 
ing rate T2 due to, e.g., inhomogcneous broadening of A. 
The laser generation requires V£ max > 7/t o r 2 . It will 
be the preferred decay mode of the inverted population 
of spin states if ujq — \JW 2 + A 2 /h coincides with one 
of the frequencies of the quantized phonon modes of the 
crystal. For a small crystal, this condition can be sat- 
isfied only at certain values of W and A. This can be 
achieved by continuously changing longitudinal or trans- 
verse magnetic field. Since the distance between these 
resonant field values is inversely proportional to the size 
of the body, the resonance condition, at a finite 7 fco , must 
be satisfied at any field for a body of size L>v t /j ko ■ 

Our results on phonon superradiance alter the existing 
theory of EPR in nanomagnets at low temperature. They 
show that one cannot properly account for the width of 
the EPR by simply computing one-phonon processes for 
isolated magnetic atoms. When the wavelength of the 
phonons, Ao = 27r/fco, is greater than the distance be- 
tween magnetic atoms, the collective processes come into 
play, greatly enhancing the rates. Observation of EPR 
lines requires fc^T < Hujq, which at T < IK translates 



into Ao < 50 nm. Thus, in the kelvin temperature range, 
our results on superradiance are applicable to nanomag- 
nets of size below 50 nm and to thin films of thickness less 
than 50 nm. Larger samples should exhibit superradiance 
at lower temperatures. The direct study of phonon su- 
perradiance requires inverted population of spin states. 
In conventional paramagnets this can be achieved in a 
pulsed magnetic field. In microscrystals of molecular 
nanomagnets the inverted population of spin states is 
easy to achieve even in a slowly varying field 4]. Note 
that in crystals of small size the phonon superradiance 
always wins over the electromagnetic superradiance [j| . 

The existence of the laser phonon mode maybe consis- 
tent with the experimental evidence [US EH that the 
low-temperature spin-phonon rates in molecular nano- 
magnets are few orders of magnitude higher than the 
computed one-molecule spin-phonon rates [ToL ITlj . In 
molecular magnets the two-level system is formed every 
time the external magnetic field, B, produces a resonance 
between the states described by the magnetic quantum 
numbers m and to'. For the (to, to') resonance, S in all 
our formulas must be replaced by \m — m!\/2, while W = 
\m — m'\gfi B B. Consider, e.g., a typical field-sweep ex- 
periment in Fes, when the field crosses the (10, —8) reso- 
nance and is 0.1T past the resonance. Taking A ~ 10 _4 K 
(l2|. M ~ 5 x 10~ 21 g, and v t ~ 10 5 cm/s, one obtains: 
luo ~ 10 n s-\ Tt ~ 10-h- 1 , but r L , max ~ 10 4 s-\ 
which would provide the to = 10 to' = —8 spin- 
phonon relaxation in less than 1-ms time. This may 
explain spin-phonon avalanches observed in sufficiently 
large crystals of molecular magnets 0|. The laser 
rate Tl = V^ o max /jk ^ as a distinctive ojq dependence 
that should be seen in experiment. At low temperature 
7 fe cx ujq due to the scattering of phonons by the in- 
homogeneities of the crystal 0, while V£ o max cx uj . 
Consequently, Tl cc ujq 3 , which is opposite to the ujq de- 
pendence of the one- molecule spin-phonon rate, Ti ocwq. 

This work has been supported by the NSF Grant No. 
EIA-0310517. 
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